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We have calculated the crystal electric field of TmTe (T > Tq) and have obtained that the ground 
state of a Tm 4/ hole is the Fy doublet in agreement with Mossbauer experiments. We study the 
quadrupole interactions arising from quantum transitions of 4/ holes of Tm. An effective attraction 
is found at the L point of the Brillouin zone, (Jl. Assuming that the quadrupolar condensation 
involves a single arm of we show that there are two variants for quadrupole ordering which are 
described by the space groups C2/c and C2/m. The Landau free energy is derived in mean-field 
theory. The phase transition is of second order. The corresponding quadrupole order parameters 
are combinations of T2g and Eg components. The obtained domain structure is in agreement with 
observations from neutron diffraction studies for TmTe. Calculated lattice distortions are found to 
be different for the two variants of quadrupole ordering. We suggest to measure lattice displacements 
in order to discriminate between those two structures. 



I. INTRODUCTION 

The monochalcogenides of rare-earth elements crystal- 
lizing in the rocksalt structure, exhibit interesting electri- 
cal, optical and magnetic propertiesJJ In TmTe thulium 
is in a divalent state with one 4/-hole (4/-'^'^) localized 
at each Tm site. At room temperature TmTe is a semi- 
conductor with the electron band energy gap 0.25 eV 
(from photoemission data)rand an / — > c? band gap 0.35 
eV (from absorption data)El. Magnetic structure deter- 
mination leads tcua type II antiferromagnet below the 
Neel temperatureQ T/v = 0.43 K. Under pressure of 
GPa TmTe undergoes a semiconductor-metal transitio: 
where a ferromagnetic order occuraJ with a Curie tem- 
perature Tc = 14 K for 2.7 GPa. At pressure of 5.7 GPa 
TmTe undergoes a structural phase transition to a tetra&. 
onal structurecnZl with a further decrease of resistivity.cl 
A renewed interest to TmTe has emerged from the re- 
cent unexpected discovery of an antiferroquadrupolar 
(AFQ) phase transition below Tq ~ 1.8 K at atmospheric 
pressure.u The traasition was first revealed by specific 
heat measurementsEl and then confirmed by elastic con- 
stant anomaliegO imd by neutron diffraction in an exter- 
nal magnetic field.lij In the quadrupolar ordered state by 
the application of a uniform magnetic field superlattice 
neutron Bragg peaks, corresponding to the wave vector 
q = (27r/a)(l/2, 1/2, 1/2), appear in addition to the nor- 
mal ferromagnetic intensity superimposed on the nuclear 
peakal3 (see for a review Ref. ^ and for recent develop- 
ment Refs. [l^ , [l3| ). Thus TmTe has been found to belong 
to an increasing family of rare-earthalil (TmCd, TmAu2, 
CeBg, DyB2C2 and so on) which exhibit quadrupolar or- 
dering. 

An important insight in the nature of the AFQ order- 
ing of TmTe can be gained on the basis of a classical 
description of multipoles presented in Ref s. [l5| , p^ . How- 
ever, the theoretical treatment in Refs. 15,16| is based on 



a phenomenological Hamiltonian where the quadrupole- 
quadrupole interactions are taken to be isotropic in the 
subspaces of Fa {Eg ) and F5 (T2g) irreducible representa- 
tions of the cubic point group Oh- In addition, a coupling 
between these representations is neglected. Such a treat- 
ment has shortcomings and does not allow to determine 
the crystal structure. Indeed, since the quadrupolar in- 
teraction forces are of short range, it is necessary to take 
into account the structure of the face centered cubic (fee) 
lattice and the full matrix character of the quadrupolar 
interactions. Furthermore, the coupling of T2g and Eg 
representations is important and affects the quadrupolar 
phase transition. It is the purpose of the present pa- 
per to investigate the quadrupolar phase transition and 
to determine the crystal structure (space group) of the 
quadrupolar ordered phase. Hence it is necessary that 
the symmetries of quadrupolar interactions on a fee lat- 
tice are carefully taken into account. Multipolar inter- 
actions of Eg and T2g modes on a fee lattice have been 
studied earlier in connections with the problern-of orien- 
tational order in a molecular crystal, solid C6oJi3 

In constructing a microscopic theory of quadrupolar 
ordering in T mTe .we will use the concepts we have devel- 
oped reeentlyt§ll3 for the description of the 7 — a tran- 
sition in Ce. We will take advantage of the fact that 
Tm and Ce are mirror elements in the series of lan- 
thanides. We recall that Ce has one single electron in 
the 4/ shell while Tm has one hole (electron configura- 
tion 4/^^). While Ce is a metal, TmTe is an insulator for 
T near Tq. Hence we will neglect the polarization of con- 
duction electrons (which contributes to the quadrupolar 
interaetionEj in Ce) in our treatment of TmTe. Because 
of charge transfer from Tm to Te, only a small fraction 
of conduction electron density is left around the Tm++ 
ion. 

In zero magnetic field the magnetic SMSceptibility 
of TmTe shows a small anomaly at TqB Thus the 
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quadrupolar ordering is not directly connected^ with the 
magnetic ordering in TmTe which occurs at stiU lower 
temperature T/y. Therefore in this paper we limit our- 
selves to the charge degrees of freedom leaving the mag- 
netic properties for future considerations. 

The content of the present paper is the following. In 
Section II wc calculate the crystal field which refers to the 
4/ hole in TmTe. The result is compared with experi- 
mental data. Next (Sect. Ill) we study the quadrupo- 
lar interactions and the resulting phase transition. We 
suggest condensation schemes leading to the monoclinic 
space groups C2/c or C2/m. In order to discriminate be- 
tween these space groups, we calculate the accompanying 
lattice distortions (Sect. IV) and suggest synchrotron ra- 
diation experiments. Finally, the results of our work are 
summarized and commented in the Conclusions (V). 



II. CRYSTAL FIELD 

The issue of the crystal electric field (CEF) of TmTe is 
controversial. A Schottky specific heat peak around 5 K 
has been found-indicating that the total CEF splitting 
is about 15 kBj From thermal expansion data at low 
temperature (2-16 K) the following sequence r8(0 K) — 
r7(10 K) - r6(16 K) has been proposed.El Recent in- 
elastic neutroa. scattering experiment supports the Fg 
ground statej£3 However, from Mossbauer spectroscopy 
another scheme, r7(0 K) - r8(12 K) - r6(19.6 K), has 
been deduced.o The latter sequence was alset,obtained 
by detailed ultrasonic velocity measurements O Here we 
would like to remark that some of these results are not 
direct and depend on conditions, methods and models 
used to fit experimental data. 

As a starting point of our derivation of the CEF we 
recall that here we consider the effect on the hole in the 
4/ shell of Tm+"'". In comparison with the case of one 
electron (as in Ce), the position of the corresponding en- 
ergy levels is reversed. For ions of lanthanides in solids 
the spin-orbit coupling 



(2.1) 



dominates crystal field effects. The potential Vso has 
spherical symmetry. Here L and S are the orbital and the 
spin angular momentum operators, J = L + S is the to- 
tal angular momentum, with eigenvalues of J^: J = 7/2 
and J = 5/2. In Eq. (2.1) ( is the spin-orbit couphng 
constant. In case of a hole, the lower energy level corre- 
sponds to the state -D7/2 (degeneracy 8) and the higher 
level to Dpj2 (degeneracy 6). The experimental energy 
separations Ago = — (7/2)C is given by 1.24 eV, corre- 
sponding to C = —4112 K. In presence of the cubic crystal 
field VcF, the hole experiences the potential 
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Since Vc f has cubic symmetry, the degeneracies p 
spin-orbit levels are lifted according to the schema^ 



(2.2) 
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TABLE I. Angular-momentum-decomposed electronic 
charges Qf' and total charges inside Tm and Te MT 
spheres and in the interstitial region (LAPW calculations, see 
Ref. 31 for details and definitions). 
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where Fg, Fg are irreducible representations of the 
double cubic group OJ^. |— ■ 

Usually for CEF calculation£3 the lowest J multiplet 
is used and the CEF 4/ Hamiltonian is expressed in-terms 
of Stevens equivalent operators Jx, Jy and Jz.EScJ A 
thorough discussion on microscopic-prigin of the crystal 
field effects is given by NewmanE3'C3 Here we will follow 
the method of Refs. |l^Jl^ which allows us to calculate 
CEF using results of electron band structure calculations 
as a starting point. Our approach is not restricted to the 
J multiplet with lowest energy. Since we may neglect 
conduction electrons (see Sect. I), the problem of the 
CEF of one 4/ hole then becomes similar to the CEF 
of one 4/ electron considered in Ref. |l^ with the only 
difference that in the following we take into account the 
radial dependence of the 4/ hole in TmTe. Our goal is 
then to diagonalize the potential and to determine 
the associated 14 eigenvalues e^. 

In order to proceed, we first need to know the charge 
distribution in the TmTe unit call. The previous elec- 
tron band structure calculationca of TmTe in the local 
density approximation (LDA) does not give such infor- 
mation. We then have performed our linear augmented 
plane wave (LAPW) calculaticms using LDA and the 
muffin-tin (MT) approximation.!^ The MT radii 2.9 and 
3.1036 a. u. were chosen for Tm and Te, respectively. The 
MT potential and density of Tm and Te have been ob- 
tained self-consistently using a LAPW basis of ~ 170 
plane wave functions on a 20-point mesh of the irre- 
ducible domain of the fee Brillouin zone. During the 
calculations the / electrons of Tm were treated as core 
states with an occupation number of 13 but were allowed 
to adjust self-consistently to the conduction electron den- 
sity. We did tabulate the radial dependence 7?./ of the 
J — 7 /2 electronic 4/ states of Tm on a set of 70 points 
with 0.111 <r< 2.971 a.u. Notice that 7?./ is obtained as 
an output of the electronic band structure calculation of 
TmTe and thereby deviates from the 4/ electronic den- 
sity of a Tm atom. The results of the charge density 
calculations are quoted in Table I. Using the MT poten- 
tial, we have calculated the spin-orbit splitting. The D5/2 
states of the 4/ hole of Tm were found to be separated 
from the /2 states by A(,q = 1.21 eV, which is close 
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to the experimental valu e of Ref. In the following we 
consider expression ( ^.l|) of Vso with = —4112 K. 

In constructing the crystal field operator Vcf, we re- 
strict ourselves to the effect of the six nearest Te neigh- 
bors in octahedral position around a Tm+"'" ion at site 
n. (Later we will discuss the limitation of the present 
approach.) Since the crystal field has cubic symmetry, 
we write it in form of a multipole expansion, following 
the method of Ref. O, Appendix A: 



TABLE II. Calculated CEF spectrum of a 4/ hole in the 
disordered phase (T > Tq), Ae=14389.4 K. 
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Here , are the hole charge density operators, Ai = 
{I — 4, Aig) refers to the lowest (non-trivial) multipole 
{I = 4) of symmetry Aig (unit representation of Oh), 
while A2 = (/ = 6, ^ig) corresponds to the next with 
I = Q. Explicitly, we have (p = 1,2): 



where the basis states \i) are 4/ states, and where 



cl,(*j) = (^|5t4lj>, (2-6a) 
cU^j)^{Ast',\j). (2.6b) 

5'aj = S^^" and = Sq^" being the site symmetry 
adapted functions^ (SAF). The function sf^" is given 
exphcitly by Eq. (A.8) of Ref. |l8|, the function S^'" is 
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The superscript F stands for quantum transitions be- 
tween 4/ states of the hole. Taking as basis states the 
functions of the irreducible representations Tiu, and 
of Oh, we obtain the diagonal coefficients {A2u) ~ 
-0.23505, cl^(T2u) = -0.03917, c^^(rij = 0.11752 
for Ai and cl^{A2u) = 0.20118, cl^{T2u) = -0.15088, 
^{Tiu) = 0.08382 for A2. Notice that there is no term 

Although in 
S'jJl, the cor- 



with ^ = 8 in the expansion series (2.4 
principle one can consider the functio 
responding matrix elements (*|S';'^^g |j) vanish for 4/ or 
bitals. The coefficients b{ and B^ are given by 
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Here eh refers to the charge of the hole, Qeff to the 
effective charge of Te and * is given by the radial 
average of the Ap-th multipole of the Coulomb potential: 



dr TZj (r) va^ o{n, n']r,r'), 
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where TZf is the radial functioiO and 
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\R{n) ~- R'{n')\' 
(2.9b) 



In the last integrals we have R{n) = X{n) + r{n) where 
X{n) is the lattice site position and f{n) the hole (elec- 
tron) radius vector, r— {r,ft), il — (8, (p). 

We take Qeff = Table I. In Eqs. (|2!9a|,b) n 

refers to a Tm lattice site, while n' - to any of its six Te 
neighbors. Diagonalization of Vcpin) leads to the crystal 
field term scheme without spin-orbit coupling: 
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^B{cl{T) + Bicl{T), 
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T2u and Ti„. 



Our calculations yield 



Bi=35.73 K, 5^=1.631 K and e{A2u) = -8.07 K, 
e(T2u) = -1.646 K, e(Ti„) = 4.336 K. Including next 
the spin-orbit coupling and proceeding as in appendix A 
of ReL [1^, we obtain the CEF term spectrum of the 
4/ hole, quoted in Table II. 

If one goes beyond the present nearest neighbor ap- 
proximation, the calculation becomes much more com- 
plex. Considering the charge contributions from the 
twelve neighboring Tm sites and Q*"* — — 1.588|e| from 
the interstitial region, we find that the crystal field split- 
tings between energy levels increase by roughly an or- 
der of magnitude in comparison with the values of Ta- 
ble II, e(r8, 1) - eiVr, 1)=62.9 K, eiTe) - e(r7, 1)=102.5 
K. The increase is mai nly c aused by a larger positive 
value of Qeff in Eqs. ( ^.84 b) when the homogeneous 
electron charge distribution in the interstitial region is 
taken into account.lEl Notice, however, that the sequence 
Fy — Fg — Fe with F7 as a ground state is conserved. 
The other remote shells of neighbors have been found to 
produce little changes on the calculation since the corre- 
sponding integrals vary likeEl \R(ji) — R{n')\~^ for S^^^ 

and \R{n) - R{n')\-'^ for S^^\ The magnitude of the 
CEF splittings in the latter approach is reduced by sev- 
eral effects: screening due to polarization of conduction 
electrons,Ej which are still present on each Tm site (sec 
Table I), and an inhomogeneous charge distribution in 
the interstitial volume. The first effect is due to the cou- 
pling of 4/ localized electrons with 5d conduotipn elec- 
tron states which was discussed by Newman.tSI In our 



3 



model the interactions with the conduction electrons are 
given by Eq. (4.10) of Ref. |9[ However, such-calcula- 
tion would require a self-consistent procedure£3 for 4/ 
and conduction electrons included in an electron band 
structure scheme and is beyond the scope of the present 
work w hijch. is focused on quadrupolar ordering. Some 
authorsBEZI have emphasized the role of covalent mixing 
which involves one-electron matrix elements between the 
4/ electrons with 5d electrons on other atoms. Notice 
that such effects imply nonoijtliogonal basis functions on 
neighboring sites. In LAPWEil and linear muffin-tin or- 
bital (LMT0)l2a methods of electron band structure cal- 
culations there is no overlap between MT basis functions 
belonging to different sites, while 4/ electronic wave func- 
tions are confined inside MT spheres. The covalent effect 
which is associated with anisotropic character of chem- 
ical bonding then may be described by inhomogeneity 
of electron charge distribution in the interstitial region 
and inside MT spheres. Since there is a large portion 
of electronic charge situated in the interstices (Q™*) it 
is likely that the effect of the charge inhomogeneity is 
appreciable. 

In summary we here obtain the sequence Fy — Fg — Fg 
with Fy as ground state, inp- agreement with results 
from Mossbauer saectroscopylHj and ultrasonic veloc- 
ity measurements. O Although the present calculation of 
CEP in incomplete our estimations show that the se- 
quence is likely to be conserved in a full treatment of the 
problem. The experimental identification of the ground 
state as Fg is in contradiction with the present calcula- 
tions. The discrepancy can arise due to strong fluctua- 
tions of quadrupole density which occurs at temperature 
T 2 — 3 K and affect the experimental results. Indeed, 
the presence of quadrupolar fluctuationsJaave been found 
at 3 K in Mossbauer studies of TmTe.Ej It is necessary 
to analyze both the experimental conditions and methods 
of obtaining the CEF levels and of identifying the states. 
Further investigations and full calculations of CEF are 
needed in order to clarify the issue. On the other hand, 
we want to stress that theoretically a quadrupole order 
can occtUp-even if the ground state of the 4/ hole is Fy 
(or F6)c3eJ while Fg remains an excited state. 



III. PHASE TRANSITION 

Here we will discuss the antiferroquadrupolar phase 
transition in TmTe. In fact this phase transition is a 
structural one, with the concomitant symmetry change 
cubic^monoclinic. In the following we will continue to 
exploit the duality electron-hole between Ce and TmTe. 
In Ref. |l^ we have shown that the Coulomb interaction 
operator between 4/ electrons (holes) on a fee lattice is 
obtained as a double multipole expansion 
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Here the expansion coefficients uaa' are given by 



vaa' {n — n) ^ J drr j dr' 

x7^^(r)7^^(/)^'AA'(n,n'; r,r'), (3.2a) 

where 
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S'a(7i) S'a(7i') 



\R{n) - R'{n')\ 

(3.2b) 



For details on the radial average in Eq. ( 3.2a ), see Ref. |19|. 
Here Sa^uX, n = (0(?i), are site symmetry adapted 

functionsEHI, A stands for (1,t), where I accounts for the 
angular dependence of the multipolar expansion and r — 
(F,fc), F denoting the irreducible representation of the 
site point group and k labeling the rows of F. In our case, 
I = 2 (quadrupoles) and F stands for the representations 
T2g {k — 1 ~3) and Eg (fc = 1, 2) of the cubic site group 
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Eqs. (2.16) and (2.15) of Ref. |l 

stands for the multipolar density 
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with 
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Introducing Fourier transforms 
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p^;(g) = -^5]e^«~-^(")pl(n), 
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VAA'{q) = Y^e'^-^^^^VAA'{h), 
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where q is the wave vector, we get for the quadrupole- 
quadrupole interaction 



t^QQ - J E E Pii'l)^^AA' (q) pL (q). 



(3.6) 
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The 5x5 matrix vaa' is given by the expressions 
(Al, A6) and (A7) of Ref. |l^. This matrix has neg- 
ative eigenvalues at some points of the Brillouin zone 
(BZ). The largest negative eigenvalues have been found 
at the X and L points of the BZ. Since the superstruc- 
ture reflectipns have been found by neutron-diffraction 
experimentalj on TmTe at the L point of the BZ, ql ~ 
(27r/a)(l/2, 1/2, 1/2), we limit in the following our con- 
siderations to the L point. 

There are four arms of the star *qL which we label 
by ql = (1/2,1/2,1/2), ql = (-1/2,1/2,1/2), = 
(1/2,-1/2,1/2) and q^ = (-1/2,-1/2,1/2), in units 
(27r/a). At q~ ql, i — 1 ~ 4:, the eigenvalue spectrum of 
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the quadrupole matrix v\/^i is the same and for simphc- 
ity we consider the arm q^. Notice that at the L point 
there is a coupUng between components of T2g and Eg 
symmetry (see (A7) of Ref. jlj) . We write 



,ET 







(3.7) 



where stands for the 2x2 zero matrix, v'^'^ describes 
the 3x3 matrix between components of symmetry, 



f3 (3 
/3 P 
/3 /3 



(3.8) 



and stands for the 3x2 T2g — Eg couphng matrix, 



A ly 
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^1 
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(^re-jt rpj^g elemen ts /3 , A, /i and v are 



and 

obtained by integrals of the type ( 3.2a ,b). 

Diagonalizing the matrix v'^'^ we obtain the eigenval- 
ues —8/3, 4/3, 4/3 and eigenvectors 



S%^{ql) = -^(SkMl) + SLi'il) + SlMD), (3.10a 



^{s'rJqD-sljqD), 



(3.10b 



S%{ql) = -^iShMl) + Sljql) - 2Sl.Jql)). (3.10c 

In the basis S%^iql), S%{qt), S%^{ql), S%^{qt), 
^E (ql) matrix v{ql) becomes block-diagonal. 



2/3 
-/3 V2i^ 
V2iy 
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We find its eigenvalues, of which X\ = — 2(— /3 + 
y/P^ + 8iy^) and A^ = 2(/3 -f ./(3^T^) are double 
degenerate while = —8/3 is non degenerate. From 
numerical calculations we obtain /3 = —33.54 K, /i = 
-29.05 K, A = 14.53 K, ly 



-224.4 K, \l 



90.3 K and A^ 



= -25.16 K and Xi = 



268.3 K. These results 
are in agreement with the symmetry relations 2iy — 
2A = —fi, that hold for quadrupole-quadrupole Coulomb 
interactions on a fee lattice. The lowest eigenvalue X\ 
has the eigenvectors 

S^'\ql) - -0.5972(^1,^^ - S^J + 0.5356 (3.12a) 
5(2) (g2) = +0.3448(5i +Sl)- 0.6895 Sl 



2g 



+0.5356 Si 



(3.12b) 





FIG. 1. Schematic pictures of two variants of quadrupole 
order parameter (quadrupole density) ; m is the mirror plane, 
C2 is the r otatio n axis, a) 5'(^' , Eg. (|3.13a|), which leads to 
C2/ c Eq. ( ^.15a| ); b) S^^', Eq. ( |3.13b[ ), which leads to C2/m, 
Eq. (|3.15b[). 



where we omit the arguments q^ on the right hand sides. 
In addition, we consider the corresponding functions in 
real space, 

+0.5556 Sl^iQ), (3.13a) 
S^'H^) = +0.3448(4,^ (f^) + S^^m 

-0.6895 S^^^in) + 0.5356 Ss^iQ). (3.13b) 

These two functions are shown in Fig. la,b. We inves- 
tigate their transformational properties in detail in Ap- 
pendix A. 

The quadrupolar densities which correspond to the 
functions S"*^"), a = 1,2, are given by the expression 



(3.14a) 



with 



dn{i\h)S'^"Hn)(h\j). (3.14b) 



(Compare with expressions (3.3), (p^).) 

The functions S'^") belong to the two dimensional small 
representation Eg of the httle group 2ni {D^d) of *qL {f^ 
representation in Kovalev's notationO). The irreducible 
representation of the space group FmSm comprises eight 
such functions, with two functions from four arms of *(fi, 
that is, 5(i)(gl), Si'Hql); 5W(g|), S^'Kql); S^'Hql), 
5(2)(g3) and S'-^^ql), S^^Hql). In principle, there are 
many possibilities for condensations-schemes at *qL in- 
volving one, two, three or four armsEj'Ej Experimentally, 
reflections associated withr-all four components of the star 
*qL were clearly observedll3 and had different intensities 
even in small applied magnetic fields. On this basis it 
was concluded in Ref. |l^ that each arm of *(Jl is asso- 
ciated with a domain. We then limit our consideration 
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FIG. 2. Monoclinic unit cell in respect to cubic system of 
axes. Black and white circles refer to Tm sites where in case 
of C2/c structure (pf ) (in case of C2/m {P2)) is taken with 
the sign + for black and — for white circles. 

to the case where a single arm, say , is involved in the 
symmetry lowering which occurs due to the quadrupolar 
ordering. In such case the following two condensation 
schemes are possibleE3: 



Fm3m : L+ [{pi (ql)) = VNpi] C2/c {Z = 2), 



(3.15a) 

Fm3m : L+ [(pf {ql)) = ^/Np2] C2/m {Z ^ 2). 

(3.15b) 

Here Paio) stands for the Fourier transform of p^{n), (...) 
denotes a thermal average and pa are the order parameter 
amplitudes. Correspondingly, in real space we obtain 



{p^{n)) = p^ cos{ql ■ X{n)), a = 1,2. 



(3.16) 



Both structures are monoclinic, with the mirror plane 
[110], see also Figs. la,b. As follows from Eq. ( |3.16| ) 
{pf{n)) and (/3|"(fi)) change from +pi to —pi and from 
+P2 to —p2, respectively, along the [110] direction. The 
resulting pattern as well as the monoclinic unit cell are 
shown in Fig. 2. We speak of an antiferro-quadrupolar 
order. At ql ther e are still three variants of condensa- 
tions of the type (3.15a) and three of the type (3.15b). 
For example, we consider the condensation in C2/c, 
Eq. ( 3.15a| ). The three variants involve condensations 
of quadrupolar functions S"^^^ which are obtained from 
5^^) through rotations by the angle 0, 27r/3 and 47r/3 
about the cubic axis [111]. Notice that this functions 
can be expressed in terms of a linear combinations of 
5(1)^ gi"^) since they form a basis of the little group 3m 
(Dsd) of ql- Otherwise, the consideration is the same as 
for S'^^\ S*^^^. The corresponding monoclinic unit cells 
are obtained from that in Fig. 2 through the same ro- 
tations by the angle 0, 27r/3 and Att/3. For example, 
for the rotation by 27r/3 the basal {xy) plane (see Fig. 
2) transforms to (yz) while the monoclinic mirror plane 



[110] becomes [Oil]. These three variants correspond to 
the so-called "S domains" whicjii,have been observed in 
neutron diffraction experimentst2l for a given arm oi*qL. 
However, the total number of domains for the condensa- 
tion to C2/c is twelve. Th e sam e holds for the second 
condensation scheme, Eq. ( ^.15b ), to the C2/m mono- 
clinic structure. We conclude that on the basis of data 
from neutron diffractionllj on TmTe it is not possible to 
deduce which of the two ordered structures actually oc- 
curs in TmTe. Both structures are monoclinic and-lead 
to the domain structure observed in experiment .Eij On 
the other hand, from our theoretical analysis of coupling 
matrices we cannot rule out one structure in favor of the 
other. In the next section we show that the two types 
of quadrupolar ordering can be distinguished by lattice 
displacements whic h acco mpany the transition. The con- 
densation scheme ( 3.15a ) corresponds to that given by 
Eq. (3.10a) of Ref. [40|, where two complex basis func- 
tions of the f ^ irreducible representation are used. 

In the following we study the thermodynamics of the 
quadrupolar phase transitions. Taking into account the 
first condensation scheme, Eq. (3.15a), we obtain for the 



quadrupole coupling at a site ni for one sublattice {{ni\) 



(3.17) 



where pi is the order parameter amplitude and where 
Pi{ni) is the quadrupolar density operator. The mean 
field Hamiltonian reads 



ff^>^^'(rM) = C/c5{,("i) + K,'(^i) 



(3.18) 



where Vj/ , Eq. (2.2), describes the crystal field and the 



spin-orbit coupling. Starting with _ff and using meth- 
ods which have been developed for molecular crystalscil, 
we obtain the following approximate expression for the 
Landau free energy: 



T/N ^Ta/N + Apl + Bp\. 



(3.19) 



Here J^o is the free energy in the disordered phase 

jc-g/iv^ -TlnZo, (3.20) 

where 



IT 



(3.21) 



is the sum of states and Si are the crys tal fi eld energy 
levels. The expansion coefficients in Eq. ( 3.19 ) are 



B = 



24[a;(2)]2 



3- 



r.(4) 



,(2)12 



(3.22a) 
(3.22b) 



where 
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TABLE III. Calculated parameters of the Landau free en- 
ergy expansion, see text for details. 





a;*-*) B/T 


k 


0.0219 


0.000768 121.5 


0.188 






(3.23a) 






(3.23b) 



ijhl 

The calculation of a;'^) and x*^^^ requires the knowledge 
of the crystal field. Using the values of from Table II 
and the corresponding calculated eigenvectors, we obtain 
the results quoted in Table III. Si nce th ere is no third or- 
der cubic invariant in expression (3.19) and since B > 0, 
the phase transition is of second order, with the transi- 
tion temperature given by 

T,^x^^^\\l\ (3.24) 

and the order parameter amplitude given at T < Tc by 



where 



k = 



12(x(2))3 



3(a;(2)) 



(3.25) 



(3.26) 



With Xl = -224.4 K we find = 4.9. This value is more 
than twice the experimental temperature Tq = 1.8 K. We 
ascribe the origin of the discrepancy to the screening ef- 
fect of conduction electrons from Tm and Te sites and 
from the interstitial region. The question may arise why 
in case of cerium the polarization of conducticm. electrons 
leads to an increase of transition temperatureEfl, while in 
case of TmTe it has the opposite effect. We recall that in 
Ce quadrupoles constructed from conduction electrons 
are in close contact and their ordering greatly reduces 
the repulsion between conduction electrons. In TmTe 4/ 
holes and conduction electrons around Tm sites are at 
larger distances and the polarization of conduction elec- 
trons merely reduces the resulting effective quadrupolar 
value. 

Finally we mention that with Eqs. ( 3.17 )-( 3.24 ) we 



readily obtain the corresp onding expressions for the sec- 
ond condensation scheme ( 3.15b ) by replacing the index 
1 in pi and ci(ij) by the index 2. The numerical values 
of x^'^\ x'^^\ B and Tc remain the same and therefore no 
distinction between C2 / c and C2 / m can be made at this 
point. 



IV. LATTICE DISTORTIONS 

The quadrupolar ordering and symmetry lowering is 
accompanied by a distortion of the cubic lattice. Such 



effects are known to occur in molecular crystals, (see for 
a review Ref. ^2]) and our present treatmentll3't£l was in- 
spired. \w the theory of orientational order in molecular 

sohds.BB 

We consider the Tm atoms located on a non rigid fee 
lattice and denote the lattice displacement of the Tm 
nucleus at site n by u{n) . For the 4/ hole coordinates we 
have 



R{n) = X{n) + r{n) + u{n), 



(4.1) 



where X{fi) stands for equilibri um nuclear position. We 
expand the intersite potential (3.1) in terms of atomic 
lattice displacements. The first order correction to the 
potential reads 

nn' 

x^(")(n)^(")(n') K(n)-u,(n')], (4.2) 

where 

(4^3, 



dX,{n) \R{n)-R'{ri')\ 



Here the index a (a = 1, 2) corresponds to the two v ari- 
ants of antiferro- quadrupole ordering, Eqs. ( 3.15a ,b). 
We recall that n — n' stands for X{n — n') = X{k). We 
take the average over the radial dependence of the 4/ 
hole: 

t;^(")(K)= / drr^ j dr' r'^ n}{r)n}{r') vI'^''\k; r,r'). 



One has the symmetry relation 



(4.4) 



(4.5) 



on the fee lattice. In the following k labels the twelve 
nearest neighboring Tm sites around a Tm ion taken as 
origin. Proceeding as in Ref. ^ we rewrite the expression 
(4.2) as an operator in the space of the 4/ hole: 



nrt' y 



xPa (") Pa (n) [ui,{ri) - M,.(n )] , 



(4.6) 



where (n) is defined by expressions ( |3.14a ,b). Trans- 
forming to Fourier space we find 

UqQT = ^EE^^'^"^(^'P>«(-^^- 9)Pa{p)Mk), (4.7) 



where 
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X cos[(p + ^) • sin[^^-^M] 



(4.8) 



In order to obtain the free energy contribution from 
Uggrp, we take the long wavelength limit fc — > and 

retain only linear terms in k. We then c onsider p near 
ql and apply the condensation schemes (3.15a,b). Fi- 
nally we replace the displacements u^{k) by their instan- 
taneous thermal expectation values {ui,{k)). After some 
algebra we obtain 



F, 



(a) 



N 



QQT 



X J2 l^k^^il^ +ikyAl;:) +ik.AiZ^]Mk)), 



(4.9) 



i^—x,y,z 

where pa is the order parameter amplitude. In ( |4.9| ) we 
have defined 

(4.10a) 

A|,"; = -aK'^'^Hl) + v'J^^ (3) + v'J"^ (4) - v'J"^ (6)], 

(4.10b) 

AiZ> = (1) + <(")(2) - <(")(4) - (5)]. 

(4.10c) 

Here the arguments k = 1 — 6 of w^*^"^ stand for X{k.), 
with Xil) = (0,1,1), ^(2) = (1,0,1), X{3) = (1,1,0), 
X(4) = (0,-1,1), X{5) = (-1,0,1), X(6) = (-1,1,0) 
in units a/2, where a is the cubic lattice constant. Intro- 
ducing the homogeneous strains 



lim ik^{uy{k)) — \/mNe^u, '^^x,y,z, 



(4.11) 



we obtain 



A(a)(r 



(4.12) 



(From the symmetry of the order parameters it follows 



that Ai"^ 



Ayy' and A^z = Ayz' .) It is convenient to 



work in the system of axes which reflects the monoclinic 
symmetry. Fig. 2. We therefore consider the coordinate 
system x'y'z', where x' axis corresponds to [110], y' - to 
[110] and z' to [001] directions of the cubic system. Notice 
that the new axes are obtained by the clockwise rotation 
about the z axis by 7r/4. Since A^^ (fJ-jV = x,y,z) is a 
tensor of the second rank, we write Fqqt [e, p] in the new 
coordinate system as 



F^^^^[e',p\/N^pl 

XX XX ^ '^^ yy y 



At^e' 



■^■'^ yz t yz\ 



(4.13) 



TABLE IV. Calculated parameters A'^i' and ho- 
mogeneous strains e^^, for the quadrupolar ordering 
Fm3m — > C2/c; pi is the order parameter amplitude. 



p V 


XX 


yy 


zz 


yz 




155.7 K 


94.2 K 


311.1 K 


243.6 K 


[^M./(P1?] X 


-2.08 


+1.95 


-4.66 


-16.0 



In the transformed coordinate system the elastic term of 
the free energy reads 

FTTieViVcN)^ \c,{e'l^ + e'ly) + ^cne'L 

+ 2c44e';„ (4.14) 

where we have introduced the notations C5 — (cn + 
ci2)/2 -I- C44, C6 = (cii -I- ci2)/2 - C44 and cn, C12, C44 
are the cubic elastic constants. Minimizing Fqqt + Ftt 
with respect to the strains e'^^ for a given configuration 
with a fixed expectation value pa we obtain: 

^'xx = -^^f|^[Al"^(ciiC5 - c?2) - A'W(cnC6 - cL) 
-2A'(^)ci2C44], (4.15a) 

^'yy = - c?2) - A'i";(cnc6 - c^) 

-2A'i^)ci2C44], (4.15b) 

^'zz = -^[A'i^)(cn + C12) - ci2(At; -f A'f-^)], 



yz — 



Pi 



2C44V^C 



(4.15c) 
(4.15d) 



where A = —20^2 + cii(cii -I- C12). The shear distortion 
e' yz implies that in the monoclinic phases the angle be- 
tween the axes y'^ and attached to the crystal deviates 
from 7r/2 by a « '^^'yz- 

We now present numerical results for the quadrupole 
order in the C2/c structure. In the monoclinic system of 
axes {x'y'z') the calculated values are quoted in Table IV. 
(The shear angle a « —32 • lO^^'p^.) For calculations of 
e'^^ we took the elastic constants cn = 10285, C12 = 3969 
and C44 = 1188 in units K/A'^ from Ref. ^ Returning 
now to the original cubic system of axes (x, y, z) we find 



-0.07- lO^Vl, e 



e^y = +2.02 • 10-Vf, 



= -4.66- 10-Vf, 
= e,y, = -11.3-10-V?- 
(4.16) 



For the condensation scheme ( 3.15b ) to the C2/m 
quadrupole structure the calculated values are quoted in 
Table V (the shear angle a = -6.84 • lO'^^pl). In the 
initial cubic system of axes we find 



8 



TABLE V. Calculated parameters A'JjJ and homogeneous 
strains e^^, for the quadrupolar ordering FmSm C2/m; p2 
is the order parameter amplitude. 





XX 


yy 


zz 


yz 




9.2 K 


488.9 K 


63.0 K 


52.3 K 




+ 12.61 


-18.87 


+1.46 


-3.42 



Starting from the quadrupole-quadrupole interactions on 
a deformable lattice, we have derived the couplings of the 
quadrupoles with the atomic lattice displacements. We 
have calculated the corresponding lattice distortions and 
suggest experiments which can be decisive in determining 
which quadrupole order is realized in TmTe. 



exx = eyy = -3.13 • IQ-^i, e,, = +1.46 • IQ-^pl 
e^y = -15.74 • lO^Vi, = = -2.42 • 10-^1- 

(4.17) 

We conclude that the two possibilities of quadrupole 
order lead to completely different displacements in the 
monoclinic phase. 
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APPENDIX A: 



V. CONCLUSIONS 

We present a microscopic model of quadrupole order 
in TmTe. For a 4/ hole above Tq we have obtained 
the sequence r7 — Fg — Fe of the crystal electric field 
(CEF) energy spectrum with Fy as ground state whichia 
in agreement with results from Mossbauer spectroscopyE^ 
and ultrasonic velocity measurements. lj The splitting of 
CEF levels is found to be small if only contributions from 
six Te nearest neighbors of a Tm site are taken into ac- 
count. 

We have considered quadrupolar interactions between 
4/ holes located on Tni. sites. On the basis of neutron 
diffraction experimental^ indicating that a single arm of 
*qL is responsible for the quadrupole structure, we have 
studied the quadrupole interactions at the L point of the 
BZ. We have found that the quadrupole coupling between 
4/ holes becomes attractive at the L point thus driving 
a structural phase transition with concomitant lowering 
of the crystal symmetry. Starting with the mean-field 
Hamiltonian we have derived the Landau free energy, cal- 
culated the transition temperature and found Tc = 4.9 K. 
The overestimation of the transition temperature is as- 
cribed to a screening effect from conduction electrons 
which has not been considered in the present work. The 
structure of TmTe below Tq is monoclinic (Fig. 2) but 
there are still two possibilities for the quadrupole or- 
der parameter. These quadrupole order parameters are 
expressed in real space in terms of T2g and Eg compo- 
nents and visualized in Fig. 1. The condensation of pf , 
Eq. (3.15a), leads to t he C2 /c structure while the con- 



densation of p2 , Eq. ( 3.15b ), leads to the C2/m. Al 



though both structures are monoclinic their symmetries 
are different. Both of them result in the domain variants 
which have been observed experimentally. We conclude 
that on the basis of the present experimental data and 
our theoretical studies it is impossible to determine un- 
ambiguously the actual quadrupole order in TmTe. We 
have shown that such discrimination could be done in 
respect to lattice distortions which develop below Tq. 



Here we investigate transforma tional properties of 
functions S'^^^Q) and S'^^^Q), Eqs. ( |3.134 b) (see Fig. 1). 
Below we omit the argument of the functions. We recall 
that the functions S^^^, S^^^ and Sj^^^ are proportional 
to the Cartesian components yz, zx and xy, respectively. 
Therefore, for the reflection m through the plane [110] 



T2,> 

2 

2g' 



while for the rotation C2 by tt about the axis [110] 



C2(S't2, 

Therefore, 



^T2o) - ('^T2„ ^T2g)^ 



2g' ^ J-23 

T2g) = +i^T2. 



l-2g 



(Al) 
(A2) 

(A3) 
(A4) 



while 



(A5) 
(A6) 



(A7) 
(A8) 
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